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I .  INTRODUCTION 


Blodynamic  computer-based  models  for  the  prediction  of  human  body 
response  to  mechanical  stress  have  become  extremely  useful  and  cost-effective 
research  and  developmental  tools,  especially  as  alternatives  to  direct 
experimentation  with  humans  and  animals.  These  models  attempt  to  simulate  or 
predict  the  forces  and  motions  experienced  by  a  body  in  high-acceleration 
events  such  as  impacts  or  from  sudden  forces  such  as  wind  shear.  In 
particular,  the  Air  Force  is  interested  in  the  reactions  of  aircrew  personnel 
to  such  forces  typically  encountered  in  various  phases  of  flight  operations, 
including  emergency  ejections  from  high-speed  aircraft.  Such  a  hazardous 
environment  is  well  suited  to  computer  modeling,  and  with  proper  execution, 
considerable  insight  into  body  motion  and  stresses  developed  in  the  body  can 
be  gained. 

The  Modelling  and  Analysis  Branch  of  the  Biodynamics  &  Bioengineering 
Division  of  the  Air  Force  Aerospace  Medical  Research  Laboratory  (AFAMRL)  has 
been  using  a  human  body  modelling  computer  program  known  as  the  Articulated 
Total  Body  (ATB)  Model  for  several  years.  The  model  is  based  on  rigid-body 
dynamics  using  Euler  equations  of  motion  with  Lagrange-type  constraints  (Fleck 
et.al.  197*4).  The  specific  configuration  uses  15  body  segments  (head,  neck, 
upper  torso,  center  torso,  upper  arms,  lower  arms,  upper  legs,  lower  legs,  and 
feet)  and  1*4  joints  between  the  segments  (Fleck  and  Butler,  1975).  Although 
it  was  originally  developed  by  the  Calspan  Corporation  for  the  study  of 
human-body  and  anthropometric-dummy  dynamics  during  automobile  crashes  for  the 
United  States  Department  of  Transportation  (Fleck  et.al.  197*4;  Fleck,  1975), 
the  ATB  Model  was  sufficiently  general  to  allow  simulation  of  whole-body 
articulated  motion  resulting  from  various  impacts  or  abrupt  accelerations 
applied  to  the  body.  Furthermore,  modifications  involving  special  Joint 


forces,  aerodynamic  forces  and  a  complex  harness  system  were  added  to 
accommodate  specific  Air  Force  applications  (Fleck  and  Butler,  1975). 

The  ATB  Model  initially  reflected  human  body  structure,  mass  distribution 
and  tissue  material  properties  for  passive  responses.  An  early  effort  to 
improve  the  ATB  Model  in  regards  to  active  responses  resulted  in  the 
development  of  a  lumped  three  parameter  viscoelastic  muscle  model  superimposed 
on  the  advance  restraint  system.  (Freivalds,  1984;  Freivalds  and  Kaleps, 

1983;  Freivalds  and  Kaleps,  1984).  However,  the  early  efforts  were 
constrained  by  the  low  number  (five)  of  harness  systems  provided  in  the  ATB 
Model,  limiting  simulations  to  simple  joint  motions  or  very  crude  whole  body 
motion.  Also,  complex  neuromuscular  functions  such  as  motor  unit  recruitment 
patterns,  time  varying  effects,  etc.  were  not  included.  Thus,  further 
development  of  the  neuromuscular  system  was  needed  to  better  simulate  active 
human  responses  to  high-g  forces. 

II.  OBJECTIVES 

The  objective  of  this  project  was  to  further  define  and  formulate 

methodologies  for  implementing  active  muscle  responses  into  the  present  ATB64«-A 

' ./  fic<\y)  fV.ath  wa  'he*’-* 

.Model. .Two  considerations  were  involved;  (1)  basic  muscle  phenomena  such  as 

C  H  J  iTA  (  w  ?  'i  V 

motor  units,  recruitment  patterns,  and  fatigue  were  to  be  included  and  (2) 
particular  emRhasis  was  to  be  placed  on  muscles  acting  in  the  torso  and  neck 
region  which  affect  flexion,  extension  and  lateral  motion  of  the  trunk  in  a 
seated  posture. 

The  objective  was  approached  in  four-phased  approach.  In  Phase  I,  the 
basic  muscle  model  developed  during  the  early  efforts  (Freivalds,  1984),  was 
re-examined  and  redefined.  In  Phase  II,  advanced  features  of  types  of  motor 
unit,  motor  unit  recruitment,  force  buildup  and  endurance  times  were  developed 
and  included  into  the  ATB  Model.  In  Phase  III,  the  ATB  Model  was  modified  to 


allow  for  up  to  50  muscles  and  a  representative  musculature  for  the  entire 
body  was  developed.  In  Phase  IV,  various  simulations  were  performed  inorder 
to  validate  the  modelling  efforts. 

III.  BACKGROUND 

A.  Skeletal  Muscle: 

Skeletal  muscles  usually  originate  on  the  skeleton,  span  one  or  more 
joints  and  insert  into  a  part  of  the  skeleton  again.  Each  muscle  is  enclosed 
in  a  connective  tissue  sheath  called  the  epimysium  and  is  held  in  its  correct 
position  in  the  body  by  layers  of  fascia.  The  muscle  is  attached  to  the  bones 
via  tendons,  while  the  interior  is  compartmentalized  into  longitudinal 
sections  called  the  fasciculi,  each  containing  many  individual  muscle  fibers. 
The  fibers  are  enveloped  by  a  connective  tissue  called  the  endomysium,  which 
transmits  the  force  of  the  muscle  contraction  from  individual  fibers  to  the 
tendons  (Fung,  1981). 

The  muscle  fibers  do  not  always  run  parallel  to  the  force  transmitting 
tendons,  as  they  do  in  fusiform  muscles.  They  can  be  arranged  in  unipennate, 
bipennate  or  multipennate  form,  thus  altering  the  force  transmittinig 
characteristics  (Fig.  1). 

The  muscle  fiber,  the  basic  structural  unit,  with  a  diameter  of  10-60  y 
and  length  from  several  millimeters  to  several  centimeters,  can  be  subdivided 
further  into  myofibrils  of  1y  diameter.  These  myofibrils  comprise  the 
hexagonal  array  of  protein  filaments  that  are  directly  reponsible  for 
the  contractile  process  and  give  rise,  with  appropriate  stains  to  the  peculiar 
striations  that  are  characteristic  of  skeletal  muscle  (Figure  2).  A  repeating 
unit  known  as  th?  sarcomere  is  defined  by  the  vertical  z-disk.  Two  types  of 
protein  filaments  are  distinguishable  in  each  sarcomere,  thin  ones  about  5nm 


UNIPENNATE  BIPENNATE  MULTI  PENNATE 

Fig.  1  Schematic  representation  of  skeletal 
muscle  fibre  arrangement. 


z-disk 


Fig.  2  Molecular  substructure  of  mammalian 
skeletal  muscle. 


(50A)  In  diameter  and  thicker  ones  about  12nm  ( 1 20 A)  across.  The  thin 
filaments  contain  actin,  globular  molecules  in  a  triple  helix,  while  the  thick 
filaments  contain  myosin,  long  molecules  with  globular  heads.  The  thin 
filaments  are  each  attached  at  one  end  to  a  z-disk  and  are  free  at  the  other 
to  interlace  with  the  thick  filaments.  The  A- band  is  the  region  of  overlap 
between  thick  and  thin  filaments,  the  I-band  contains  solely  the  thin 
filaments,  while  the  H-band  is  the  middle  region  of  the  A-band  into  which  the 
actin  filaments  have  not  penetrated  (Fung,  1981). 

The  actual  contractile  process  takes  place  at  the  junctions  between  the 
myosin  and  actin  in  a  process  known  as  the  sliding  filament  theory  first 
presented  by  H.  E.  Huxley  (1953)*  The  myosin  molecules  consist  of  a  long  tail 
piece  and  a  "head".  The  tails  lie  parallel  in  a  bundle  to  form  the  core  of 
the  thick  filament  while  the  heads  project  laterally  from  the  filament  in 
pairs,  rotated  with  respect  to  its  neighbors  to  form  a  spiral  pattern  along 
the  filament.  These  heads  seem  to  be  able  to  nod;  they  lie  close  to  their 
parent  filament  in  relaxation,  but  stick  out  to  actin  filaments  when  excited. 
Thus,  during  muscle  contraction  the  muscle  fiber  shortens  as  the  filaments 
slide  over  each  other,  forming,  breaking  and  reforming  chemical  bonds  between 
the  myosin  heads  and  the  globular  actin  molecules. 

B.  Previous  Neuromuscular  Modelling  Efforts: 

Previous  modelling  efforts  (Freivalds,  1984;  Freivalds  and  Kaleps,  1983; 
Freivalds  and  Kaleps,  1984)  produced  the  lumped  model  of  skeletal  muscle  shown 
in  Figure  3.  Structures  which  lie  in  parallel  to  the  force  producing 
sarcomeres:  the  sarcolemna  (sheath)  of  the  individual  fiber  and  the  various 
outer  connective  sheaths  (fascia,  endomysia,  perimysia)  are  represented  by  the 
parallel  elastic  element  (PE).  Practically  all  the  tension  observed  when 
stretching  the  resting  muscle  will  result  from  this  element.  Because  the 


12 


muscle  fiber  is  60f  water,  an  appropriate  damping  element  (DE)  is  also 
Included  parallel  to  the  contractile  elements  (CE).  The  contractile  element 
represents  the  purely  contractile  protein  molecules.  In  series  with  the 
contractile  element  is  a  bridge  element  (BE)  representing  the  elastic  elements 
within  the  cross  bridges  and  the  z-disks.  The  parallel  elastic  element  for 
the  sarcomere  (PS)  does  not  contain  a  damping  component,  since  the  sarcolemna 
attached  to  the  z  disks  does  not  allow  appreciable  movement.  The  tendinous 
parts  of  the  muscle  fiber  are  located  near  the  origin  and  insertion  of  the 
fiber  and  thus  are  depicted  by  a  series  elastic  element  (SE).  Any  mass  of  the 
sarcomeres  is  disregarded,  especially  when  compared  to  the  much  larger 
external  mass  that  the  muscle  contraction  must  move  (Hatze,  1981). 

The  functional  form  of  a  harness  within  the  advanced  restraint  system 
(Butler  and  Fleck,  1980)  was  utilized  in  representing  the  muscle 

F(e,e' )  -  Ft  ( e)  ♦  F2(e)F3(e')  ♦  F^e')  (1) 

where  e-strain,  e'-strain  rate  and  F-total  force.  Obviously,  from  the 
functional  form  of  Equation  1 ,  only  three  elements  can  be  modelled  adequately. 
Furthermore,  all  three  of  these  have  to  be  in  parallel,  since  forces  in 
viscoelastic  theory  add  directly  only  in  parallel.  Series  elements  would  have 
required  complex  integro-differential  equations  which  could  not  have 
replicated  with  out  extensive  changes  to  the  model.  However,  certain 
assumtions  allow  a  fairly  easy  reduction  of  the  lumped  model  to  a  simpler  form 
fitted  by  Equation  1.  SE  and  BE  can  be  considered  to  be  very  stiff  springs 
and  eliminated  completely.  This  contention  is  supported  by  Bawa  et.al.  (1976) 
who  found  KgE"3724  N/m  to  be  much  large  than  Kpg-1000  N/m.  Kgg  can  be 
considered  to  be  in  a  similar  range  with  Kse*  Eliminating  SE  and  BE  results 


In  a  model  with  four  parallel  elastic  elements.  These  can  be  combined  Into 
one  parallel  elastic  element,  yielding  the  final  simplified  model  In  Figure  4. 

The  total  force  developed  by  the  simplified  model  of  Figure  4  can  now  be 
expressed  as: 

F  -  (fpE  ♦  fCE  +  fpE^MAX  <2) 

where  F^ax  is  the  maximum  Isometric  tension  of  the  muscle.  Equation  2 
corresponds  very  nicely  with  Equation  1  with: 

fPE  FMAX  “  F1(e) 

fCE  FMAX  “  F2^  ^3) 

fDE  FMAX  “  F4(e) 


where 


F2(e)  -  fa(e) 
F3(e)  -  fy(n) 


(4) 


For  concentric  or  shortening  contractions  fpg  is  zero  (PE  producing  force  only 
under  stretch)  while  fpE  produces  a  force  opposite  in  sign  to  the  contractile 
force.  For  eccentric  or  lengthening  contractions,  fpg  and  fj)E  all  act  in  the 
same  direction  as  f(jE* 

Mathematical  representations  for  each  element  were  determined  as  follows. 
For  the  parallel  elastic  element,  extensive  tests  on  the  tensile  properties  of 
resting  human  sartorlus  muscle  carried  out  by  Yamada  (1970)  indicate  an 


exponential  force-strain  function: 


fPE  -  .ooi6296(e7*66l6e-i )  (5) 

where  fpE  is  the  force  developed  by  the  PE  normalized  with  respect  to  maximum 
isometric  tension  in  the  muscle  and  e  is  the  strain: 


where  l  is  the  instananeous  muscle  length  and  lQ  is  the  resting  length.  This 
force-strain  curve  is  shown  in  Figure  5. 

The  velocity  dependence  of  the  damping  element  (DE)  can  be  expressed 
similarly  to  the  form  used  for  a  simple  mechanical  dashpot: 


fpE  *  .00588  e 


where  fpE  is  the  normalized  force  and  e  is  the  muscle  strain  rate.  This 
curve  is  shown  in  Figure  6. 

The  contractile  element  is  the  only  active  component  in  the  model.  Its 
behavior  is  extremely  complex  and  depends  nonlinearly  on  its  length, 
contractive  history,  velocity  of  movement,  the  degree  of  stimulation  and  its 
temperature.  However,  for  practical  purposes,  only  the  basic  functions  were 
considered:  the  length-force  relationship  and  the  force-velocity 
relationship. 

The  length-force  relationship  is  determined  by  the  number  of  active  cross 
links  or  filamentary  overlap  and  can  be  adequately  expressed  from  the  data  of 
Gordon  et.al.  (1966)  by  the  functional  suggested  by  Hatze  (1981,  p.  42): 


fj,(e)  -  .32  ♦  .71e-1-112e  sin(3.722(e+.344)) 


This  function  is  shown  in  Figure  7. 


The  force-velocity  relationship  is  determined  by  the  rate  of  breaking  and 
reforming  the  cross  bridges  with  higher  rates  producing  less  effective  bonds. 
To  account  for  the  whole  range  of  negative  velocities  (shortening  or 
concentric  contractions)  as  well  as  positive  velocities  (lengthening  or 
eccentric  contractions)  Hatze  (1981,  p.  45-46)  has  defined  the  following 
expression: 

fy(n)  -  .1^33  {.1073  ♦  e"1,409  sinh  (3.2n  ♦  1 .6)  }-1  -.005[2-e6e]  (9) 

where  fv(n)  is  the  normalized  force  due  to  the  force-velocity  relationship 
as  defined  by  the  first  term  and  reduced  by  internal  resistance  as  defined  by 
the  second  term.  However,  since  the  coefficient  of  the  second  term  is  smaller 
by  a  factor  of  30  than  the  first  term,  it  can  be  disregarded  for  present 
purposes,  n  represents  the  normalized  contractile  element 
velocity: 

n  -  (to) 

with  e^AX  being  the  maximum  shortening  velocity  of  the  contractile  element. 
Equation  6  represented  by  Figure  8. 

IV.  PHASE  I  -  REDEFINITION  OF  THE  BASIC  MUSCLE  MODEL 

A.  Passive  Viscoelastic  Elements: 

Previous  efforts  to  define  a  neuromuscular  response  with  the  ATB  Model 
utilized  simulations  of  elbow-flexion  (Frelvalds,  1984;  Freivalds  and  Kaleps, 
1983).  These  indicated  several  conclusions.  The  passive  strain  function  as 
represented  by  fpg  produced  an  insignificant  increment  to  the  total  muscle 
force  even  with  15 i  strain  at  complete  extension  of  the  forearm.  This  was 


B.  Active  State  Function: 


The  action  of  the  contractile  element  is  more  complex  than  initially 
modelled.  Not  only  is  the  force  production  dependent  on  the  length  of  the 
muscle  fiber  and  the  velocity  of  contraction  but  also  on  the  active  state  q  of 
the  muscle  fiber.  This  active  state  q  is  defined  as  the  relative  amount  of  Ca 
bound  to  troponin  (inhibitor  molecule  of  actin).  If  the  maximum  number  of 
potential  interactive  sites  on  the  thin  filament  are  exposed  by  the  action  of 
Ca++,  then  q-1 ;  while  in  a  resting  state  q-q0.  Thus  the  isometric  tension 
developed  by  a  muscle  fiber  at  a  given  length  lq  the  CE  is  directly 
proportional  to  q  (Hatze,  1981,  p.  33)* 

Define  Y  to  be  the  difference  between  the  real  free  Ca++  concentration  Yf 
and  the  free  Ca++  concentration  Y0  in  the  resting  fiber.  However,  for 
practical  purposes  since  Y0  <<  Yf,  we  have  Y-Yf.  Let  p  -  dq/dY  denote  the 
Ca++  concentration  rate  of  change  of  the  active  state  q.  The  process  of 
binding  Ca++  ions  to  the  troponin  sites  i3  hypothesized  by  Hatze  (1981)  and 
supported  by  the  experimental  studies  of  Ebashi  and  Endo  (1968)  to  be  a 
function  of  the  length  l  of  the  CE  and  of  the  difference  between  the  maximum 
and  present  value  of  q  and  controlled  by  a  negative  feedback  loop  as 
follows. 

dp/dY  -  Pi2(e)  (1-q)-2p2Pi (e)p  (8) 

where  e  -  (l-t0)/t0  is  the  strain  and 

dq/dY  -  p  (9) 

Solving  the  differential  system  of  Equation  8  and  Equation  9  with  initial 
conditions: 


p(0)  -  0 


above  50 %  strain.  Reductions  in  force  levels  due  to  the  active  length-force 
relationship  were  larger  amounting  to  approximately  a  1 6 %  reduction  of  the 
maximum  voluntary  contraction  (MVC)  at  full  flexion  and  a  10%  reduction  at 
full  extension.  The  force  velocity  curve  showed  as  much  as  30%  changes  in 
force  production  levels. 

The  force  produced  by  the  passive  viscous  damping  element  amounted  to  as 
much  as  a  30%  increase  above  MVC.  This  seemed  to  be  an  excesesively  large 
effect  for  a  passive  response  and  consequently  the  damping  coefficient  was 
re-evaluated.  Glantz  (1974)  derived  a  damping  coefficient  of  5(g/mm2)/(mra/sec) 
using  the  data  of  isolated  cat  papillary  (cardiac)  muscle  of  length  8.5mm  from 
Parmley  et.al  (1970).  This  value  was  normalized  to  a  more  useful  unitless 
value  of  .00588  expressed  as  fractional  force  per  velocity  of  muscle 
lengths/second.  Force  was  normalized  by  the  maximum  muscle  strength  of  100  N 
per  cm2  of  muscle  cross-sectional  area.  This  seems  to  be  the  upper  limit  of 
inherent  muscle  strength  as  determined  by  a  wide  range  of  researchers  using 
different  techniques  (Fick,  1910;  Ikai  and  Fukunga,  1968;  Hatze,  1981). 

Velocity  was  normalized  by  the  muscle  specimen  length  of  8.5mm  such  that  the 
fDE  calculation  could  be  used  regardless  of  muscle  type  or  length.  This 
unknown  factor  in  the  previous  deviation  (Freivalds,  1984)  was  probably  the 
main  cause  of  the  unusually  large  values  of  damping  forces.  A  final 
assumption  necessary  to  complete  the  calculation  of  the  damping  coefficients 
is  that  the  coefficient  for  skeletal  muscle  would  be  similar  in  value  to  that 
of  cardiac  muscle.  Since  coefficients  for  skeletal  muscle  have  not  been 
calculated,  no  other  alternative  is  available. 


r.*-*-  -*  - 


q(0)  -  q0  -  .005 


One  obtains  a  normalized  solution: 


,  v.  ,  (’-0o>  f  "2»1 

<><C*T>  ■  I  V 


m2P1(e)Y  (e)Y 


where 


,  2  1\1/2  S  1 

p2  i  (p2  -  D  p2  >  i 


“1  ,2  K2  -  'K2 


Substituting  experimentally  found  values  (Hatze,  1981): 


p2  -  2.34  x  10^  (^77^")  -.*»*»  <  e  <  -8 


P2  ■  1.05 


One  obtains: 


q(„Y)  -  1  -  (1-qo)(2.U  »<c)Y( ,,,-2.096  X  107h(e)Y(t)) 


where 


e+.l  ' 


However,  a  simple  computational  approximation  is  provided  by  Hatze  (1981,  p. 
40)  for  most  mammualian  muscles: 


q(e,Y) 


qo+p2(e)Y2 

1+p2(e)Y2 


where 


*  .  •  ■  •  •  .  •  •  *  -"•,*••• 

* •  •  •  •  * 


,  v  ^  1 .9(e+1 ) 

p(e)  -  66,200  -  ■ 


(16) 


The  function  Y,  the  free  Ca**  ion  concentration  can  be  represented  as  a 
function  of  time  t  and  the  stimulation  rate  v  by  a  trend  function  which 
represents  the  average  behavior  of  Y  in  successive  time  intervals,  and  which 
approaches  a  maximum  value  asymptotically  and  has  the  rate  of  increase 
proportional  to  the  stimulation  rate  (Hatze,  1981,  p.39). 

Y  -  m(Cv-Y)  Y(0)  -  Y0  (17) 

where  m  and  c  are  constants  (C  -  1  .373  x  10**  (Hatze,  1981)  and  m  to  be 
determined  later)  and  v  is  the  relative  stimulation  rate  defined  by 

0  <  v  -  -  <1  (18) 

“  T  " 

where  t-1  and  t_1  denote  the  stimulation  rate  and  maximum 
stimulation  rate  respectively. 

The  rate  of  stimulation  of  motor  units  during  voluntary  contraction  has 
been  very  controversial.  Several  studies  have  found  a  fairly  constant 
discharge  frequency  over  a  wide  range  of  tension  for  individual  motor  units 
(Bigland  and  Lippold,  1954;  Clamann,  1970),  while  others  maintain  that  an 
increase  in  muscle  tension  is  achieved  in  part  by  an  increase  in  the 
stimulation  rate  and  that  this  may  be  important  in  achieving  precision  and 
smoothness  of  contraction  (Marsden  et.al  1971;  Person  and  Kudina,  1972; 
Milner-Brown  et.al  1973).  However,  even  for  those  studies  who  found  a  rise  in 
discharge  frequency  with  tension,  the  frequency  at  the  start  of  the  discharge 
for  rapid  contractions  was  much  higher  and  closer  to  the  maximum  stimulation 
rate  (Tanji  and  Kato,  1973).  Thus,  it  is  fairly  reasonable  to  assume  a 
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constant  stimulation  rate  and,  therefore,  a  constant  relative  stimulation 


rate. 

Now  solving  Equation  17  with  v  constant  yields: 

Y  -  (Y0-Cv)  e"fflt  +  Cv 

With  Y0  <<  Cv  (Y0  -  1x10^9;  Hatze,  1 98 1 ) ,  Equation  19  reduces  to: 

Y  -  Cv(1-e"fflt) 


(19) 


(20) 


Substituting  Equation  20  and  Equation  16  into  Equation  15  yields  (Figure  9): 


q(t) 


.005  ♦  82.63  v2  (1-e~mt) 
1  ♦  82.63  v2  (1-e~mt) 


(21) 


Consequently,  fcE  can  be  redefined  by  using  the  relative  force  fq  developed  by 
the  active  state  function  q: 

f  -  f  (t)  f  (c)  f  (n)  (22) 

oh  q  p  v 

IV.  PHASE  II  -  ADVANCED  DEVELOPMENT 

A.  Organization  of  Fibers  into  Motor  Units 

The  population  of  motor  units  can  be  subdivided  into  two  distinct 
populations  based  on  their  contractile  and  histochemical  properties:  Type  I 
(slow  twitch)  motor  units  and  Type  II  (fast  twitch)  motor  units  (Close,  1972). 
Type  I  motor  units  have  slower  contraction  times,  tend  to  be  more  aerobic  and 
less  fatigable  and  are  recruited  at  lower  tension  levels.  Type  II  motor  units 
have  faster  contraction  times,  tend  to  be  more  anaerobic  and  more  fatigable 
and  are  recruited  at  higher  tension  levels  (Close,  1972;  Milner  -  Brown, 


Fig.  9  RELATIVE  FORCE  DEVELOPED 
BY  ACTIVE  STATE  FUNCTION 
VS.  TIME 
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et.al.,  1973).  Thus,  the  total  muscle  force  output  should  be  the  sum  of  the 
force  output  from  Nj  of  the  Type  I  motor  units  and  Nji  of  Type  II  motor 
units: 


(fPE  +  fCE, 


+  f 


+  f  )F 

CEn  DE'  MAX 


(23) 


Furthermore,  the  same  total  population  can  be  subdivided  into  two-dynamlcally 
different  populations:  the  N  population  of  active  motor  units  and  the  N-N 
population  of  inactive  or  resting  motor  units,  where  N  is  the  total  number 
of  motor  units  in  the  muscle. 

Muscle  properties  dependent  on  fiber  type  will  be  developed  in  the  next 
section  as  many  of  these  also  depend  on  the  recruitment  pattern  used. 

Orderly  Recruitment  of  Motor  Units 

It  has  been  well  established  that  motor  units  are  recruited  in  a 
sequential  order  according  to  their  sizes  (Milner  -  Brown,  et.al.  1973).  The 
cumulative  relative  cross-sectional  area  u  occupied  by  the  fibers  of  the 
recruited  units  increases  by: 

u  -  u  e°  — —  o  <  u  <  u  <  1  (24) 

o  -  ° 

where 

c  -  -in  u  , 
o 

N  is  the  number  of  stimulated  motor  units  and  N  is  the  total  number 
of  motor  units  (Hatze,  1979).  For  N  large,  Equation  24  reduces 


(25) 


a 

n 


cN/N 

e 

o 


a  (1-N/N) 
o 


l  (1-n) 
o 


where  n  Is  the  normalized  number  of  recruited  or  active  motor  units  -  N/N. 

The  Aju  of  the  relative  cross-sectional  area  u  upon  recruitment  of  the 
1th  motor  unit  is  then  defined  by: 


AjU  -  C  e°i/N  (26) 

where  C  is  a  normalization  constant  determined  by  the  requirement  that 
N 

that  £  Aiu  -  1  i.e. 

1-1 


c  .  - } -  (27) 

N  c,/N 

l  e1 
1-1 

Applying  the  ratio  of  the  smallest  to  the  largest  motor  unit  potential 
measured  in  a  muscle  to  Equation  25,  an  estimate  of  the  value  of  uQ  for  a 
given  muscle  can  be  found.  These  range  from  uo-.005  for  the  human  rectus 
femorus  muscle  to  uo«.009  for  the  human  biceps  muscle  (Hatze,  1979)  with  an 
average  value  of  u0«. 00673  resulting  in  c-5  to  be  used  for  the  present  study. 

Thus  two  very  important  properties  of  motor  unit  recruitment  dynamics 
have  been  included;  motor  units  are  normally  recruited  sequentially  from  the 
smallest  to  the  largest  and  the  size  of  the  recruited  units  grown 
exponentially.  Combining  the  two  sets  of  overlapping  population  distribution 
of  motor  units  yields  two  distinct  cases:  a)  N  <  Nj,  ie.  only  part  of  the 
Type  I  motor  units  are  stimulated  and  none  of  the  Type  II  can  be  stimulated 
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because  of  the  orderly  recruitment  pattern  and  b)  N  >  NIt  ie.,  all  Type  1 
motor  unite  are  stimulated  and  some  of  the  Type  II  motor  units  are  stimulated, 
but  none  of  the  Type  I  are  inactive.  These  conditions  can  be  expressed 
mathematically  by  adjusting  fgg  by  the  motor  unit  type  and  the  relative  area 
of  the  number  of  motor  units  recruited.  Since  the  force  developed  in  a  muscle 
is  directly  proportional  to  the  area  stimulated  (Ikai,  M.  and  Fukunga,  T. 


N. 


1968).  Thus  for  0  <  n  <  n^.  where  n^. 


a)  fCEt  ‘  f,I(t)  flI<E)  V">  Uo 


(i-n) 


(28) 


and  for  n^.  <  n  <  1 


b)  -  f  (t)  f.  (e)  f  (n)  u o0""!5 

9t  *i  v! 


CE, 


(29) 


f  -  f  (t)  f  (e)  f  (n)  Cu  (1_n)-u 0(1'nI)3 
CEII  qII  *-11  VII  °  ° 


Equation  28  and  29  can  be  substituted  directly  into  Equation  23  for  the 
overall  definition  of  muscle  force,  fpg  and  f^g  are  determined  by  the  total 
cross-sectional  area  of  the  muscle,  whether  all  motor  units  are  completely 
recruited  or  not.  One  further  adjustment  is  to  substitute 

fMAX  “  kA  (30) 

into  Equation  23.  The  maximum  isometric  tension  for  any  skeletal  muscle  is 
determined  by  the  inherent  muscle  strength  per  cross  sectional  area  k  (k  - 
100  N/cm^  as  determined  by  Pick,  1910;  Ikai  and  Fukunga,  1968;  Hatze,  1981) 
multiplied  by  the  cross-sectional  area  of  the  muscle  A  (Figure  10). 


Additional  properties  of  motor  units  can  be  obtained  from  the  data 
presented  in  Henneoan  and  Olson  (1965).  The  contraction  time  tc  of  a  motor 
unit  is  a  decreasing  function  of  the  fraction  n  of  recruited  motor  units: 

tc  -  a2~a3n  (31) 

Thus  for  Type  I  motor  units 

tcl  ~  *2^*3^  0  <  nx  <  nx  (32) 

and  for  Type  II  motor  units 

t,  -  a0  -  a-  n  nT  <  n  <  1  (33) 

CII  2II  3II  1 


The  constants  a  a-,  a  _  ,  a,  can  be  determined  from  experimental 
i  51  II  ^II 

values.  For  n-0,  the  value  of  tc  corresponds  to  the  contraction  time  of  the 
slowest  Type  I  unit  in  the  muscle,  approximatley  equal  to  .1  sec.;  for  n-n* 
the  value  of  tc  corresponds  to  the  fastest  Type  I  unit,  approximately  .045 
sec.;  and  for  n-1  (given  nj*1 )  the  value  of  tc  corresponds  to  the  fastest  Type 
II  unit,  approximately  .02  sec.  (Stephens  and  Stuart,  1975).  Substituting  and 
solving  for  the  unknown  values  yields: 


a  -  .055/n 

1 


.02  + 


.025 

(1-nj) 


.025 

O-nj) 


(34) 
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Substituting  Equation  31*  into  Equation  32,  33  yields 


.  .  .055  n 

tr  -  .1  -  — 

CI  nI 


.02  + 


'II 


.025  Hi 
1-n, 


(35) 

(36) 


Several  other  important  parameters  can  be  derived  using  Equation  31 . 

Close  (1965)  showed  that  for  mammalian  skeletal  muscle  the  maximum  normalized 
speed  of  shortening  is  related  to  the  contraction  time  of  a  muscle,  consisting 
predominantly  of  one  fiber  type,  by: 


emax 


(37) 


where  B  has  a  value  of  .297  for  human  muscle.  Using  Equation  35  and  36, 
cmax  13  found  to  be  2.97/sec.  for  slow  (n-0)  and  14.85/sec.  for  fast  (n-1 ) 
motor  units.  However,  for  present  modelling  purposes,  an  average  value  for 
cmax  wil1  be  used.  Integrating  tc  from  Equation  35  over  the  pattern  of  motor 
of  motor  unit  recruitment: 


J 

0 


yields: 


(.1  - 


^55x> 


r  ox. 

J  e  dx 


.1  (ecn-1)  - 


.055 
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ecn(cn-1  )-1 


cn  , 
e  -1 


(38) 


(39) 
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A  Similar  process  Is  used  for  t  in  Equation  36. 

CII 

The  rate  constant  in  Equation  19  also  depends  on  the  contraction  time: 

m  -  A0/tc  (40) 

where  A0  -  .372  for  human  muscle  (Hatze,  1981,  p.  62)  For  slow  motor  units 
(n-0)  m  -  3-72  and  for  fast  motor  units  (n-1)  m  -  18.5.  Again  an  average 
value  is  used  for  present  modelling  efforts. 

Time  Varying  Effects 

The  most  important  time  varying  effect  in  the  muscle  is  fatigue.  It  is 
obvious  that  people  can  maintain  their  maximum  effort  very  briefly  (5 
seconds),  whereas  they  can  maintain  a  force  of  around  a  quarter  of  their 
maximum  strength  for  an  extended  period  of  time.  Such  am  endurance  responses 
can  be  explained  by  examining  the  properties  of  individual  motor  units.  Type 
I  motor  units  tend  to  be  more  aerobic,  less  fatigable  and  are  recruited  at 
lower  tension  levels.  While  Type  II  motor  units  tend  to  be  an  aerobic,  more 
fatigable. and  are  recruited  at  higher  tension  levels,  (Stephens  and  Usherwood, 
1977).  Although  exact  fatigue  and  recovery  patterns  for  individual  motor 
units  have  not  been  identified,  the  maximum  endurance  time  cam  be  estimated 
from  experimental  studies.  The  earliest  experiments  of  Miller  (1932)  implied 
that  the  length  of  time  a  force  could  be  maintained  depended  on  the  fraction 
of  available  strength  to  be  exerted.  This  relationship  was  further  verified 
by  Rohmer t,  1960,  Kogi  and  Hakamada,  1962;  Caldwell  1963.  1964;  Monod  and 
Scherrer,  1965;  Schutz,  1972.  Only  three  studies  attempted  to  derive  and 
publish  formulas  of  this  relationship.  Monod  and  Scherrer  (1965) 
proposed: 


34 


Kogl  and  Hakamada  (1962)  suggested 


E(mln) 


5012 


(*FMAX} 


1.99 


(42) 


while  Schutz  (1972)  Indicated: 


E(mln) 


■1 .25  ♦  1 .25 


*F 


MAX 


(43) 


All  of  these  formulas  have  some  faults  that  limit  their  usefulness  in 
representing  the  empirical  data.  Equation  42  and  43  do  not  account  well  for 
the  asymptotic  relationship  of  endurance  approaching  indefinite  times  for 
force  levels  of  15-20*  MVC,  Equation  41  does  provide  the  asymptote  but 
predicts  lower  than  normal  endurance  times  for  large  force  levels.  A  separate 
formula  was  developed  for  the  current  work,  based  on  the  data  of  Rohmert 
(I960),  who,  with  over  300  subjects  tested,  had  the  largest  sample  size.  Best 
fit  was  produced  by  the  hyperbolic  relationship  (shown  in  Figure  11): 


E(sec) 


1236.5 

(*fmax'15) 


72.5 


(44) 


Once  the  endurance  time  is  exceeded,  however,  the  person's  strength  does 
not  immediately  fall  to  zero.  For  maximal  or  large  submaximal  efforts,  there 
is  still  a  gradual  decay  to  the  lower  level  of  15-20  percent  found  for 
indefinite  holds  (Petrofsky,  1982,  p.  55).  This  experimental  data  can  be 
modelled  very  easily  using  polynomial  regression: 


%F  -  98.1  -23. 9t  ♦  1.9t2 


(45) 


VI.  PHASE  III  -  MODELLING  THE  GENERAL  MUSCULATURE 

The  attachment  of  the  complete  muscle  systems  to  limb  segments.  Includes 
the  Identification  of  Joint  biomechanics,  the  measurement  of  origin  and 
Insertion  distances,  the  integration  of  agonist  and  antagonist  actions  and 
computation  of  cross-sectional  areas  for  estimation  of  total  force  production. 
Methods  to  accomplish  this  can  be  best  described  according  to  the  joints  or 
area  of  the  body  involved. 

A.  Elbow  Joint: 

Modelling  of  the  elbow  and  simulation  of  elbow  flexion  is  perhaps  the 
easiest  case  to  examine  and  will  serve  as  a  simple  examole  demonstrating  the 
validity  of  the  technique  used  for  the  more  complicated  joints.  The  model 
Includes  two  segments:  the  humerus,  and  a  combination  of  the  ulna  and  the 
radius  and  three  elbow  flexor  muscles:  biceps  brachil,  brachialls  and 
brachioradlalls.  Examination  of  the  biceps  brachil  more  closely  shows  the 
origin  of  the  long  head  to  be  at  or  beyond  the  gleno  -  humeral  joint  (He  Minn 
and  Hutchings,  1977).  The  insertion  can  be  set  at  approximately  3.5  cm  from 
the  elbow  Joint  corresponding  to  the  data  of  Wilkie  (1950).  Using  the 
cross-sectional  area  of  4.58  cm2  for  the  biceps  (Schumacher  and  Wolff,  1966) 
and  multiplying  by  the  maximum  muscle  force  of  100N/cra2  (Hatze,  1981)  yields  a 
maximum  isometric  tension  of  458  N.  (The  difference  between  anatomic  and 
physiological  areas  due  to  fiber-  orientation  are  accunted  for).  Adding  the 
force  of  690N  generated  by  the  brachialls  (6.9  cm2  x  100  n/cm2)  to  the  biceps 
yields  a  maximum  elbow  flexion  force  of  1148  N.  Such  a  value  can  be  compared 
to  the  data  of  Wilkie  (1950)  who  found  that  his  subjects  could  maintain  a 
maximum  of  195.8  N  at  the  wrist.  With  a  lever  ratio  between  muscle  insertion 
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distance  and  the  moment  ar<o  of  the  weight  of  .15,  the  maximum  elbow  flexion 
force  is  1305  N.  These  values  are  remarkably  close  considering  that  other 
minor  muscles  producing  additional  torgue  are  not  accounted  for  in  the  first 
calculation.  Similar  calculations  were  conducted  on  the  other  muscles. 

Complete  details  on  all  of  the  elbow  muscles  are  based  on  previously 
collected  data  (Chao  and  Morrey,  1978;  Youm,  et.al.  1979;  Maton,  et.al.  1980; 
Amis,  et.al.,  1980;  Hatze,  1981)  and  are  summarized  in  Table  1. 

B.  Shoulder  Joint: 

The  shoulder  is  a  .much  more  complicated  joint  consisting  of  three 
separate  Joints:  the  glenohumeral  Joint,  the  acromioclavicular  Joint  and  the 
sternoclavicular  Joint.  Correspondingly,  many  more  muscles  are  involved  to 
produce  many  different  actions.  Details  on  the  actions  of  these  muscles, 
points  of  origin  and  insertion,  along  with  cross-sectional  areas  are  taken 
from  previous  biomechanical  studies  (Schumacher  and  Wolff,  1966;  Dempster, 
1965;  DeLuca  and  Forrest,  1973;  Engln,  1980)  and  given  in  Table  2. 

C.  Hip  Joint  and  Knee  Joint: 

Although  consisting  of  only  one  joint,  the  hip  is  a  ball  and  socket  joint 
and  along  with  the  many  muscles  Involved,  undergoes  many  different  types  of 
actions.  A  further  complication  is  that  a  majority  of  these  muscles  span  both 
hip  and  knee  Joints  and  thus  hip  actions  cannot  be  uniquely  separated  form 
knee  actions.  Details  on  the  actions  of  these  muscles,  points  of  origin  and 
Insertion,  along  with  the  cross-sectional  areas  were  taken  from  earlier 
biomechanical  studies  (Merchant,  1965;  Schumacher  and  Wolff,  1966;  Seireg  and 
Arvlkar,  1973;  Jensen  and  Davy,  1975;  Crowninshleld,  et.al.  1978;  Dostal  and 
Andrews,  1981;  Smldt,  1973;  Nissan,  1980;  Wlsmanls,  1980;  Minns,  1980)  and  are 
summarized  in  Table  3> 
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D.  Trunk  and  Neck  Musculature: 


Simulation  of  the  trunk  musculature  is  a  much  more  difficult  undertaking 
than  for  the  previous  Joints.  First  of  all,  there  are  many  muscles  involved, 
close  to  20  major  ones  for  the  lower  back  and  trunk  and  equally  many  for  the 
neck  region.  Secondly,  some  of  the  muscles,  such  as  the  longus,  spinalis  and 
semispinalis,  have  many  attachment  sites  between  the  different  vertebrae. 
Thirdly,  the  lines  of  action  of  the  muscle  forces  are  not  always  in  straight 
line,  e.g.,  the  interior  and  exterior  obliques.  Fourthly,  the  vertebral 
Joints  are  complicated  by  the  ligaments  and  their  additional  force-bearing 
capabilities.  Appropriate  approximations  were  used  when  necessary. 

Details  on  insertions  and  origins  were  obtained  from  Rab  et.al.  (1977), 

Rab  (1979),  Takashima,  et.al.  (1979),  and  Williams  and  Belytschko  ( 1 98 1 ) ,  while 
cross-sectional  areas  were  used  from  Schumacher  and  Wolff  (1966)  and  Williams 
and  Belytschko  (1981).  A  summary  of  these  findings  is  given  in  Tables  4 
(neck)  and  5(trunk). 


Table  1  Specifications  on  Elbow  Musculature 


Cross- 

Sectional 


Muscle  Group 

Origin 

Insertion 

Action 

Area  (cm2) 

Biceps  brachii 

Short  head  from 
coracoid  process 
of  scapula,  long 
head  from 
supraglenoid 
tuberosity 

Radial  tuberosity 

Flexion  of 
forearm 

3.55 

Brachial is 

Lower  anterior 
surface  of 
humerus 

Coronoid  tuber¬ 
osity  of  ulna 

Flexion  of 
forearm 

H.63 

Brachiora- 

dialis 

Proximal  two- 
thirds  of 
humerus 

Styloid  process 
of  ulna 

Flexion  of 
forearm 

1.37 

Triceps 

brachii 

Long  head  from 
infraglenoid 
tunerosity  of 
scapula,  lateral 
and  medial  head 
from  posturior 
surface  of  humerus 

Olecranon 

Extension  of 
forearm 

16.38 

Anconeus 

Lateral  epicondylel 
humerous  ' 

Olecranon 

Extension  of 
forearm 

.9»1 

Pronator 

teres 

Medial  eplcondyle 
of  humerus 

i 

Middle  of  radius 

Pro nation 
of  forearm 

1.61 

Supinator 

Lateral 
epicondyle  of 
humerus 

Later ial  and 
anterior  surface 
of  radius 

Supination 
of  forearm 

1.77 

Table  2  Specifications  on  Shoulder  Musculature 


Cross- 

Sectional 


Muscle  Group 

Origin 

Insertion 

Action 

Area  (cm2] 

Deltoid 

Clavicle,  scapula 
acromion 

Deltoid  tuberosity 
of  humerus 

Abduction 
of  arm 

11  .01 

Suprasplnatus 

Supraspinous  fossa 
of  scapula 

Greater  tuberosity 
of  humerus 

Abduction 
of  arm 

3.3 

Pectoralis 

major 

Clavicle,  Stermum 

Bicipital  groove 
of  humerus 

Adduction 
of  arm 

6.8 

Latissimus 
dor  si 

Lower  Thoracic 
and  lumbar 
vertebrae 

Bicipital  groove 
of  humerus 

Adduction 
of  arm 

5.37 

Teres  major 

Inferior  angle  of 
scapula 

Bicipital  groove 
of  humerus 

Adduction 
of  arm 

H.97 

Teres  minor 

Axillary  border 
scapula 

Greater  tuberosity) 
of  humerus 

Adduction 
of  arm 

1.57 

Subscapular  is 

Subscapular  fossa 
scapula 

Lesser  tuberosity  j 
of  humerus 

Flexion 
of  arm 

9.9 

Coracobrachl 

alls 

Coracoid  process 
scapula 

Medial  border  of 
humerus 

Flexion  of 

arm 

1.52 

Infraspinatus 

Inf rasp lnous 
fossa  of  scapula 

Greater  tuberosity 
of  humerus 

Extension 
of  arm 

5.98 

Table  3  Specifications  on  Hip  and  Knee  Musculature 


Cross- 

Sectional 


Muscle  Group 

Origin 

Insertion 

Action 

Area  (cm2) 

Gluteus 

medius 

I.4 1  r  ^ 

ISaHH 

Greater  trochanter 
of  femur 

Abduction 
of  thigh 

21 .18 

Gluteus 

minimus 

Outer  surface  of 
ilium  (pelvis) 

Greater  trochanter 
of  femur 

Abduction 
of  thigh 

9.6 

Tensor  fasiae 
latae 

Anterior  part  of 
iliac  crest  of 
pelvis 

Iliotibial  tract 
of  femur 

Abduction 
flexion  of 
thigh 

2.48 

Obturatus 

internus 

Obturatus  foramen 
area  of  pelvis 

Greater  trochanter 
of  femur 

Abduction 
of  thigh 

3-91 

Adductor 

longus 

Pubis  of  pelvis 

L,inea  aspera  of 
femur 

Adduction 
flexion  of 
thigh 

5.03 

Adductor 

brevis 

Pubis  of  pelvis 

Below  lesser 
trochanter  of 
femur 

Adduction 
flexion  of 
thigh 

4.54 

Adductor 

magnus 

Ischial 
tuberosity  of 
pelvis 

Linea  aspera  of 
femur 

Adduction 
flexion  of 
thigh 

20.58 

Pectineus 

Pubic  tubercle 
of  pelvis 

Below  lesser 
trochanter  of 
femur 

Adduction 
flexion  of 
thigh 

2.47 

Quadratus 
f  eraor  is 

Ischial 
tuberosity  of 
pelvis 

Quadrate  tubercle 
of  femur 

Adduction 
of  thigh 

2.91 

Obturatus 

externus 

Obturator  foramen 
area  of  pelvis 

Trochanteric 
fossa  of  femur 

Adduction 
of  thigh 

4.95 

Gluteus 

maximus 

Iliac  crest  of 
sacrum 

Iliotibial  tract 
of  femur 

Extension 
of  thigh 

29.42 

Semimem¬ 

branosus 

Ischial  tuberosity 
of  pelvis 

Upper  part  of 
tibia 

■ 

12.97 

Semi tend i- 
nosus 

Ischial  tuberosity 
of  pelvis 

Medial  condyle 
tibia 

III 

4.33 

Quadriceps 
f emor is 


Iliopsoas 


Gastrocnemus 


Ischial 
tuberosity  of 
pelvis,  linea 
aspera  of  femur 

Lateral  condyle 
of  tibia,  head 
of  fibula 

Iliac  spine  of 
pelvis  anterion 
surface  of 
femur 

Patela 

L2-L4  vertebral 
bodies,  iliac 
fossa  of  pelvis  • 

Lesser  trochanter 
of  femur 

Medial  and 
lateral  condyles 
of  femur 

Calcaneus 

Lateral  condyle 
of  femur 

Posterior  surface 
of  tibia 

Pubic  symphysis 
of  pelvis 

Upper  medial 
surface  of 
tibia 

Iliac  Notch  of 
of  pelvis 

Upper  Medial 
surface  of 
tibia 

Longus  Capitis 


Rectus  capitis 
an ter  lor 


Rectus  capitis 
lateralis 


Rectus  capitis 
posterior  major 


Rectus  capitis 
posterior  minor 


Obllquus 

capitis 

superior 


Longlsslmus 

capitis 


Spinalis 

capitis 


Transverse 
processes  of  C3. 
C4,  C5 ,  C6 

Basilar  part  of 
occipital  bone 

Flexion  of 
head 

Cl  transverse 
processes 

Front  of  foramen 
magmum  on 
occipital  bone 

Flexion  and 
notation  of 
head 

Cl  transverse 
processes 

Jugular  process 
of  occipital 
bone 

Lateral 
flexion  of 
head 

C2  spinous 
process 

Inferior  nuchal 
line  of  occipital 
bone 

Extension, 
lateral 
flexion  of 
head 

Cl  spinous 
process 

Inferior  nuchal 
line  of  occipital 
bone 

Extension, 
later! al 
flexion  of 
head 

Cl  transverse 
processes 

Inferior  nuchal 
line  of  occipital 
bone 

Extension 
and  lateral 
rotation  of 
head 

Spinous  process 
of  T1  and  C7 

Occipital  bone 
and  temporal 
bone 

Extension 
and  lateral 
flexion  of 
head 

Transverse 
processes  of 

T1 ,  C6,  C4 

Mastoid  process 
of  Temporal 
bone 

Extension 
and  Lateral 
flexion  of 
head 

Transverse  pro¬ 
cesses  of  T1  and 
C7 

Between  superior 
and  inferior 
nuchal  line  of 
occipital  bone 

Extension 
and  lateral 
flexion  of 
head 

Transverse  pro¬ 
cesses  of  T1  and 
C7 

Between  superior 
and  Inferior 
nuchal  line  of 
occipital  bone 

Extension 
and  lateral 
flexion  of 
head 

Table  4  (Continued) 


Trapezius 

Heads  of 
clavicles  and 
and  spines  of 
scapulae 

Occipital  bone 
and  thoracie 
vertebrae 

Extension 
lateral 
flexion  of 
head 

10.6 

Sternocleidomas- 

toldeus 

Head  of  sternum, 
medial  sections 
and  heads  of 
clavicles 

Occipital  bone 
and  temporal 
bone 

Lateral 
flexion  and 
flexion  of 
head 

1.6 

Levator  scapulae 

Medical  sections 
of  scapulae 

Transverse  pro¬ 
cesses  of  Cl ,  C3 

Lateral 
flexion  of 
neck 

17.75 

Longus  colli 

Anterior  side  of 
body  of  C5 
anterior  side  of 
body  of  C6 
anterior  side  of 
body  of  T1 

Anterior  side  of 
body  of  C4 
anterior  side  of 
body  of  C3 
anterior  side  of 
body  of  C4 

Flexion  of 
neck 

.75 

Scalenus 
anterior,  mediua 
and  posterior 

Medial  clavicle 

C3»  C4,  C5,  C6, 

C7  transverse 
processes 

Flexion  and 
later ial 
flexion  of 
neck 

1.75 

Splenius 

cervicis 

T1  spinous  pro¬ 
cess 

Transverse  pro¬ 
cesses  of  Cl ,  C2 

C3 

Lateral 
flexion  of 
neck 

.7 

Longissimus 

cervicis 

T1  transverse 
processes 

Transverse  pro¬ 
cesses  of  C2,  C3» 
C4 ,  C5 ,  C6 

Extension  of 
neck 

.6 

Spinalis 

cervicis 

Spinous  process 
of  T1  and  C7 

C2  spinous  process 

Extension  of 
neck 

1 .25 

Semispinalis 

cervicis 

Transverse  pro¬ 
cesses  of  T2  and 
C7 

Spinous  process 
of  C2,  C3,  C4  and 
C5 

Extension 
of  neck 

2.00 

Table  5  Specifications  on  Trunk  Musculature 


Cross- 

Sectional 


Muscle  Group 

Origin 

Insertion 

Action 

Area  (cm2) 

Illocastalls 

lumborum 

Transverse  pro¬ 
cesses  L1-L5 

Lower  six  ribs 

Extension 
and  lateral 
flexion  of 
vertebral 
column 

1.0 

Iliocostalis 

dorsl 

Lower  six  ribs 

Upper  six  ribs 

Extension 
and  lateral 
flexion  of 
vertebral 
column 

.5 

Longlsslmus 

dorsl 

Transverse  pro¬ 
cesses  L1-L5 

Transverse  pro¬ 
cesses  T1-T12 

Extension 
and  lateral 
flexion  of 
vertebral 
column 

1.0 

Spinalis  dorsl 

Spinous  processes 
L2,  LI,  T1 2,  Til 

Spinous  processes 
TU-T8 

Extension 
and  lateral 
flexion  of 
vertebral 
column 

1 .0 

Semlspinalls 

dorsl 

Transverse  pro¬ 
cesses  T7-T12 

Spinous  processes 
C6,  C7,  T1 ,  T2, 

T3 

Extension 
and  lateral 
flexion  of 
vertebral 
column 

1.0 

Multlf ldus 

Transverse 
processes  C5-T12 

.. 

Spinous  processes 
above  vertebra  of 
origin 

Extension 
and  lateral 
flexion  of 
vertebral 
column 

Intersplnales 

Spinous  pro¬ 
cesses  L5-C2 

Spinous  processes 
above  vertebra 
of  origin 

Extension  off 

vertebral 

column 

.5 

Intertrans- 
versar 11 

Connect  adjacent 

transverse  processes 

Lateral 
flexion  of 
vertebral 
column 

.25 

Table  5  (Continued) 


External  Anterior  half  of  Lower  8  ribs  Compresses  6.8! 

Oblique  lilac  crest  _  abdomen 


Internal  oblique 

Anterior  half  of 
iliac  crest 

Lower  3  ribs 
and  mid  line  of 
body 

Compresses 

abdomen 

5.68 

Rectus 

abdominus 

Pubic  symphysis 

Xyphoid  process 

Flexes 

vertebral 

column 

2.66 

Quadratus 

Lutnborum 


Iliac  crest 


Transverse  pro¬ 
cesses  L1-L4 


Flexes 

vertebral 

column 
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VII  PHASE  IV  -  SIMULATION  AND  VALIDATION 

In  Phase  IV,  the  advanced  near ornuse alar  model  was  validated  via  simalatlon 
of  haman  body  responses  to  high  G( lateral)  acceleration.  Data  obtained  ander 
similar  conditions  on  air  crew  personnel  experiencing  tests  in  the  Dynamic 
Enviornmental  Simaaltor  at  Wright-Patterson  AFB,  OH,  were  ased  for  comparison 
purposes.  Some  of  this  data  was  used  in  an  earlier  pilot  study  developing  a 
trunk  musculature  (Frelvalds,  1984). 

The  full  ATB  Model  with  15  body  segments  (head,  neck,  upper  torso,  center 
torso,  pelvis,  upper  arms,  lower  arms,  upper  legs,  lower  legs, and  feet)  and  14 
Joints  (head  Junction,  neck  junction,  waist,  L5/S1,  Joint,  two  hips,  two  knees 
two  ankles,  two  shoulders,  two  elbows)  was  utilized  to  provide  an  adequate  human 
neuromusculature.  The  muscles  described  in  Chapter  VI  and  summarized  in  Tables 
1-5  needed  to  be  added  to  the  existing  ATB  Model.  For  each  muscle  specific 
coordinates  of  origin  and  insertion  points  were  determined  form  anatomical  texts 
(Qulrlng  et.al.,  1945s  Quirlng,  1947;  Gray,  1974;  McMinn  and  Hutchings,  T 9 T7 ) . 
Cross  sectional  areas  obtained  from  Schumacher  and  Wolff  (1966)  and  Williams  and 
Belytschlso  (1981)  were  multiplied  by  the  muscular  force  constant  of  100  N/cm^ 
(Ikai  and  Fukunaga,  1968)  to  determine  a  muscle  force  scaling  factor.  These 
resulting  values  were  then  converted  to  English  units  for  use  in  the  present  ATB 
Model  mode  and  are  summarized  in  Tables  6-10. 

For  the  simulations  of  responses  to  lateral  G  forces,  the  body  segments 
were  arranged  in  the  semi -reclining  posture  maintained  by  air  crew  personnel  in 
the  cockpit.  The  lower  trunk  was  restrained  by  a  lap  belt;  any  other  restrains, 
such  as  shoulder  pasds  or  hands  placed  on  controls,  were  eliminated.  Only  the 
neck  and  trunk  musculature  on  the  right  side  of  the  body  were  activated  so  as  to 
reduce  program  complexity  and  execution  times.  A  2  Gy  lateral  force  was  applied 
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to  the  body  and  the  acceleration;  velocity  and  displacement  of  various  body 
segments  were  recorded. 

A  graphical  response  of  the  whole  body  reponse  (using  only  trunk  and  neck 
musculature)  to  the  lateral  force  over  time  is  shown  in  Figure  13.  For 
comparison  purposes,  the  response  to  a  control  case  with  no  musculature  is  given 
in  Figure  12a,  while  the  response  using  the  previous  simplified  musculature  is 
given  in  Figure  12b.  Although  the  musculature  does  not  completely  prevent  the 
lateral  deflection  of  the  body,  the  response  is  significantly  delayed  with  head 
and  neck  maintaining  the  upright  position  for  a  longer  period  of  time.  The 
result  is  better  observed  in  Figure  14,  which  shows  the  plot  of  angular 
displacement  of  the  upper  trunk  for  allthree  conditions.  At  the  end  of  256  msec 
the  angular  displacement  is  reduced  by  20°  with  the  use  of  musculature.  Based 
on  the  time  history,  the  response  with  the  musculature  lags  upto  40  msec  behind 


the  control  response. 


Table  6  ATB  Specifications  for  Elbow  Musculature  (Right  Side  Only) 


Force 

ORIGIN  (in)  INSERTION  (in)  Scaling 

Factor 


Muscle  Group 

SEG 

X 

Y 

Z 

SEG 

X 

Y 

Z 

(lbs) 

1) 

Biceps  -  Short  Head 

3 

1.57 

6.28 

-3.24 

13 

.39 

.0 

-6.56 

35.7 

Brachii  -  Long  Head 

3 

.0 

6.28 

-2.85 

13 

.39 

.0 

-6.56 

42.4 

2) 

Brachial is 

12 

.39 

.0 

-2.68 

13 

.39 

.0 

-7.15 

101.86 

3) 

Brachioradialis 

12 

.39 

.39 

3.61 

13 

.0 

.79 

.33 

30.14 

4) 

Triceps  Medial  Head 

12 

-.39 

.0 

-  .72 

13 

-.59 

.0 

-8.53 

116.9 

Brachii  Lateral  Head 

12 

-.39 

.0 

-1.5 

13 

-.59 

.0 

-8.53 

115.0 

Long  Head 

3 

.O'- 

6.88 

-1.85 

13 

-.59 

.0 

-8.53 

128.4 

5) 

Anconeus 

12 

.0 

1.18 

5.19 

13 

.39 

.0 

-6.94 

20.68 

6) 

Pronator  Teres 

12 

.39 

.0 

4.35 

13 

.0 

.71 

-3.81 

35.42 

7) 

Supinator 

12 

.0 

1.18 

5.19 

13 

.0 

.59 

-5.94 

38.94 

Table  7  ATB  Specifications  for  Shoulder  Musculature  (Right  Side  Only) 


ORIGIN  (in)  INSERTION  (In)  Scaling 

Factor 


Muscle  Group 

SEG 

X 

Y 

Z 

SEG 

X 

Y 

Z 

(lbs) 

1  ) 

Deltoid  a) 

3 

-.79 

6.49 

-3.81 

12 

.39 

.39 

.07 

121  .11 

b) 

3 

.79 

1  6.49 

-3-81 

12 

.39 

.39 

.07 

121.11 

2) 

Suprasplnatus 

3 

.0 

3.74 

-1.85 

12 

.0 

1.57 

-5.83 

72.6 

3) 

Pectoralls  a) 

3 

2.0 

.0 

-.28 

12 

.39 

.39 

-1 .9 

74.8 

Major  b) 

3 

2.0 

3.34 

-2.64 

12 

.39 

.39 

-1.9 

74.8 

4) 

Latissimus  a) 

2 

-2.0 

.0 

3.85 

12 

.0 

.39 

-2.68 

59.18 

Dorsl  b) 

3 

-2.0 

.0 

.71 

12 

.0 

.39 

-2.68 

59.18 

5) 

Teres  Major 

3 

.0- 

4.92 

2.09 

12 

.0 

-.39 

-3.08 

109.34 

6) 

Teres  Minor 

3 

.0 

5.71 

.51 

12 

.0 

.39 

-5.44 

34.54 

7) 

Subscapular  Is 

3 

.39 

4.92 

.51 

12 

.59 

.59 

-5.44 

217.8 

8) 

Coraco  Brachial Is 

3 

1.57 

7.28 

-3.23 

12 

.0  1 

-.39 

-.72 

33.44 

9) 

Infraspinatus 

3 

-.39 

4.52 

-.28 

12 

.0  j 

.59 

-5.44 

131 .56 

K 


I  .• 

P 

»V 

I".* 

F 


w".* 

tv 

tv 

C* 


Table  8  ATB  Specifications  for  Hip  and  Knee  Musculature  (Right  Side  Only) 


ORIGIN  (in) 


INSERTION  (in) 


Scaling 

Factor 


Muscle  Group 

SEd 

X 

2 

SEG 

_ xJ 

Y 

2 

(lbs) 

1) 

Gluteus  Medius 

1 

-1 .46 

.0 

-2.23 

6 

BIS 

5.88 

-8.21 

465.96 

2) 

Gluteus  Minimus 

1 

.0 

Ha 

-1.44 

6 

-1.18 

5.88 

-8.21 

211  .2 

3) 

Tensor  Fasciae  Latae 

1 

-1.57 

-2.23 

6 

-1.97 

5.49 

-5.26 

54.56 

4) 

Obturatus  Internus 

1 

.0 

1.95 

3.68 

6 

-1.18 

3.49 

-8.88 

86.02 

5) 

Adductor  Longus 

1 

.0 

.76 

3.28 

6 

-1 .1 8 

4.7 

-1.32 

110.66 

6) 

Adductor  Brevis 

1 

.0 

1.55 

3.68 

6 

-1.1  8 

4.7 

-3-09 

99.88 

7) 

Adductor  Magnus  a) 

1 

.0 

3-13 

3.68 

6 

— i .  i  a 

4.7 

-6.44 

226.38 

b) 

1 

.0-. 

1.55 

3.68 

6 

-1 .1  8 

4.7 

-8.8 

226.38 

8) 

Pectineus 

1 

.787 

1.55 

.92 

6 

-1.18 

1.0 

-5.06 

54.34 

9) 

Quadratus  Femoris 

1 

-.79 

1.95 

3.68 

6 

-1.97 

-7.62 

64.02 

10) 

Obturatus  Externus 

1 

.0 

1.95 

3.68 

6 

-1.18 

-1.97 

-8.8 

108.9 

11) 

Gluteus  Maximus  a) 

1 

.0 

-3*15 

-.65 

6 

-1.1  fi 

1.38 

-4.47 

323.62 

b) 

1 

.0 

.0 

.0 

6 

-1 .18 

1.38 

-4.47 

323.62 

12) 

Semimembranosus 

1 

-1.57 

2.49 

7 

.0 

-1.18 

-7.14 

285.34 

13) 

Semi tend inosus 

1 

-1.97 

.0 

2.49 

7 

.39 

-.79 

-4.39 

95.26 

14) 

Biceps  Femoris  a) 

1 

-1.77 

i .  i  a 

2.49 

7 

.0 

1.57 

-7.14 

129.8 

b) 

6 

-1.18 

1.18 

-.53 

7 

.0 

1.57 

-7.14 

129.8 

15) 

Quadriceps  a) 

1 

.0 

-.23 

7 

1.18 

.0 

-7.93 

308.0 

b) 

6 

.0 

1.18 

.65 

7 

1 .18 

.0 

-7.93 

924.0 

16) 

Iliopsoas  a) 

1 

.0 

-.39 

3.41 

6 

-1.10 

.79 

-6.44 

165.66 

b) 

2 

.0 

.0 

-.17 

6 

-1.18 

.79 

-6.44 

165.66 

17) 

Gastrocnemius 

6 

.0 

.0 

8.52 

8 

.0, 

.0 

-1.5 

344.52 

18) 

Popliteus 

6 

-1.18 

-1.57 

9.7 

7 

.0 

-5.57 

43.78 

19) 

Gracilis 

1 

-.39 

HE 

3.67 

7 

.34 

-.79 

-5.17 

35.86 

20) 

Sartorius 

1 

-1.57 

.0 

-1.05 

7 

.39 

-.79 

-5.17 

34.1 

%  .*% 


%  / 


// 


Table  9  ATB  Specifications  for  Neck  Musculature  (Right  Side  Only) 


ORIGIN  (in)  INSERTION  (in)  Scaling 

Factor 


Muscle  Group 

SEd 

X 

Y 

z 

SEG 

X 

Y 

Z 

(lbs) 

1) 

Longus  Capitis 

4 

-.75 

1.0 

-.54 

5 

1.54 

.39 

3.18 

16.5 

2) 

Rectus  Capitis 
Anterior 

4 

-.75 

1 .0 

-2.52 

5 

1 .26 

.39 

3.18 

5.5 

3) 

Rectus  Capitis 
Lateralis 

4 

-.75 

1.0 

-2.52 

5 

.87 

1.18 

3.18 

5.5 

4) 

Rectus  Capitis 
Posterior  Major 

4 

-.75 

.0 

-1.86 

5 

-.55 

.94 

2.83 

11.0 

5) 

Rectus  Capitis 
Posterior  Minor 

4 

-.75“ 

.0 

-2.52 

5 

-.55 

.31 

2.98 

8.47 

6) 

Obliquus  capitis 
Superior 

4 

-.75 

1 .0 

-2.52 

5 

-.55 

1.38 

2.67 

22.0 

7) 

Splenius  Capitis 

4 

-.75 

.0 

2.16 

5 

-.39 

1.57 

2.67 

26.84 

8) 

Longissimus  Capitis 

4 

-.75 

1.0 

2.16 

5 

.08 

1 .34 

2.67 

11.0 

9) 

Spinalis  Capitis 

4 

-.75 

1.0 

2.16 

5 

-.87 

.24 

2.83 

11.0 

10) 

Semispinalis  Capitis 

4 

-.75 

1.0 

2.16 

5 

-1 .42 

.79 

2.67 

52.36 

ID 

Trapezius 

3 

.0 

7.28 

-2.35 

5 

-1.42 

.39 

2.59 

233.2 

12) 

Sternocleido- 
Masto ideus 

3 

3.0 

.0 

-2.0 

5 

-.39 

1.57 

2.67 

35.2 

13) 

Levator  Scapulae 

3 

-2.0 

4.28 

-1 .85 

2 

-.31 

.45 

-1 .85 

390.5 

14) 

Longus  Colli 

3 

.75 

.0 

-2.0 

4 

.75 

.0 

-.18 

16.5 

15) 

Scalenus 

3 

3-0 

2.0 

-1.85 

4 

-.75 

1.0 

.2 

38.5 

16) 

Splenius  Servicis 

3 

-.75 

.0 

-1.26 

4 

-.75 

1 .0 

-1.86 

15.4 

17) 

Longissimus  Servicis 

3 

-.75 

1.0 

-1.26 

4 

-.75 

1.0 

-.54 

13.2 

18) 

Spinalis  Cervicis 

3 

-.75 

.0 

-2.6 

4 

-.75 

.0 

-1.86 

27.5 

19) 

SemiSpinalis 

Cervicis 

3 

-.75 

1.0 

-2.6 

4 

-.75 

.0 

-1.2 

44.0 
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Table  10  ATB  Specifications  for  Trank  Mascalatare  (Right  Side  Only) 


ORIGIN  (in)  INSERTION  (in)  Scaling 

Factor 


Muscle  Group 

SEd 

X 

Y 

Z 

SEG 

X 

Y 

Z 

(lbs) 

1) 

Iliocostal is 

Luoborum 

1 

-1.0 

1.5 

3.85 

2 

-1.0 

3.0 

.31 

22.0 

2) 

Iliocostalis  Dorsi 

2 

-1.0 

3.0 

.31 

3 

-1.0 

3.0 

.0 

11.0 

3) 

Longisslmus  Dorsi 

2 

-1  .0 

1.5 

-.19 

3 

-.75 

1.0 

.71 

22.0 

4) 

Spinalis  Dorsi 

2 

-1.0 

.0 

-4.23 

3 

-.75 

1.0 

.71 

22.0 

5) 

Semispinalis  Dorsi 

3 

-.75 

1.0 

1.89 

4 

-.75 

1.0 

1.5 

22.0 

6) 

Multifidus  a) 

2 

-1 .0 

1.5 

-2.89 

3 

-.75 

0 

4.23 

12.75 

b) 

3 

-.75 

1.0 

-2.04 

4 

-.75 

0 

2.82 

13-75 

7) 

Interspinales  a) 

2 

1.0 

.0 

-2.89 

3 

-.75 

0 

4.23 

5.5 

b) 

3 

-.75 

.0 

-2.04 

4 

-.75 

0 

2.82 

5.5 

8) 

Intertrans-  a) 

Versarii 

2 

-1.0 

1.5 

-2.89 

3 

-.75 

1.0 

4.23 

2.75 

b) 

3 

-.75 

1.0 

-2.04 

4 

-.75 

1.0 

2.82 

2.75 

9) 

External  Oblique 

1 

3.0 

3.52 

-3.08 

2 

.0 

3.52 

.31 

150.7 

10) 

Internal  Oblique 

1 

.0 

3.52 

-3.08 

2 

3.0 

3.52 

.31 

124.95 

11 ) 

Rectus  Abdominus 

1 

1.0 

.5 

1.92 

3 

3.0 

.5 

-4.23 

58.52 

12) 

Quadratus  Lumborum 

1 

.0 

3.52 

-3.08 

2 

-.75 

1 .0 

-.17  ! 

61.6 
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ANGULAR  DISPLACEMENT  (ROLL)  OF  UPPER 
TRUNK  DUE  TO  LATERAL  G-FORCES 


TIME  -  mt«c 


VIII.  CONCLUSIONS 


Further  simulations  of  various  human  body  motions  or  responses  to  external 
forces  need  to  be  conducted  in  order  to  more  adequately  validate  the  model. 
These  would  include  not  only  additional  lateral  forces,  but  also  forward/back 
(Gx)  and  up/down  (Gz)  forces  as  well  as  multidirectional  forces.  Various 
combinations  of  muscle  parameter  values  need  to  be  tested,  so  as  to  determine 
the  optimum  values  for  generating  the  most  realistic  human  response. 

Similarly,  further  examination  of  neuromuscular  factors  such  as  reflex  arcs 
may  still  lead  to  further  Improvements  In  the  model.  However,  the  simulations 
demonstrated  that  the  presently  developed  muscle  model  can  adequately 
represent  an  active  human  neuromusculature  response  to  dynamic  mechanical 
stresses  and  can  serve  as  a  cost  effective  research  and  developmental  tool. 
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